Sensitivity of microwave transition in H2O2 to variation of the electron-to-proton mass 

ratio 
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Recent observation of several microwave transitions in H2O2 from the interstellar medium 
[Bergman et al. Astron. Astrophys., 531, L8 (2011)] raised interest to this molecule as yet an- 
other sensitive probe of the tentative variation of the electron-to-proton mass ratio /x. We estimate 
sensitivity coefficients of the microwave transitions in H2O2 to /x- variation. The largest coefficient 
for 14.8 GHz transition is equal to 37, which is comparable to highest sensitivities in methanol and 
an order of magnitude higher than sensitivity of the tunneling transition in ammonia. 

PACS numbers: 06.20.Jr, 06.30.Ft, 33.20.Bx, 95.85.Bh 



Molecules with an internal motion of large ampli- 
tude recently attracted much attention as high sensi- 
tive probes of the tentative variation of the electron- 
to-proton mass ratio /x = m e /m p . In 2004 van Veld- 
hoven et al. [lj pointed out that inversion transition in 
ammonia is very sensitive to //-variation. Later, sev- 
eral groups used this transition to place stringent lim- 
its on /^-variation on a large timescale of the order 
of 10 Gyr 0-11] • At the same time ammonia spectra 
from our Galaxy were used to study possible depen- 
dence of fi on the local mater density [g, 0|- More re- 
cently the sensitivity to /x- variation was studied for other 
molecules with tunneling motion, including hydronium 
(H 3 0+) [|,IH, methanol (CH 3 OH) [IM1, and methy- 
lamine (CH3NH2) [Hj]. All these molecules are observed 
in the interstellar medium (ISM) and potentially can be 
used as probes of /x- variation. Sensitivities of the tunnel- 
ing transitions are two orders of magnitude higher com- 
pared to optical transitions in molecular hydrogen, which 
is traditionally used to study /x- variation at high redshifts 
(see [3] and references therein). Even higher sensitivi- 
ties can be found for certain mixed tunneling-rotational 
transitions. Recently methanol and methylamine were 
detected at rcdshift 2 = 0.89 First limit on /x- 

variation using methanol is reported in (Tlj . 

Peroxide molecule (H2O2) is one of the simplest 
molecules with large amplitude tunneling mode. It is 
very well studied both theoretically and experimentally. 
However, to the best of our knowledge, there is no anal- 
ysis of the sensitivity of the microwave transitions in this 
molecule to /x-variation. This is probably explained by 
the fact that H2O2 was not observed in ISM. Situation 
has changed after the first observation of four microwave 
lines of H2O2 from the molecular cloud core p Oph A 
in our Galaxy by Bergman et al. [16]. In this communi- 
cation we estimate sensitivities of the microwave mixed 
tunneling-rotational transitions to /x-variation. 

In equilibrium geometry H2O2 is not flat; the angle 27 
between two HOO planes is close to 113°. Two flat config- 
urations correspond to local maxima of potential energy; 
the potential barrier for trans configuration (27 = ir) is 
significantly lower, than for cis configuration (7 = 0): 
w 400 cm" 1 and U « 2500 cm" 1 . To a first ap- 



proximation one can neglect the tunneling through the 
higher barrier. In this model peroxide is described by 
the slightly asymmetric oblate top with inversion tunnel- 
ing mode, similar to ammonia and hydronium. A more 
accurate theory accounts for tunneling through both bar- 
riers [13] • In this case torsion motion can be described 
as hindered rotation. For sufficiently high angular quan- 
tum numbers J and Ka this internal motion strongly 
interacts with overall rotation. According to [l8[ this in- 
teraction becomes important for Ka > 7. Such levels lie 
very high and can not be observed from ISM. Below 70K 
there are only levels with quantum numbers J < 6 and 
Ka < 2. For such levels the simpler model is sufficiently 
accurate and we will use it here. 

When tunneling through both barriers is taken into ac- 
count the ground torsion state splits in four components 
designated by the quantum number t = 1 - 4. Because 
trans barrier is much lower than the cis one, the split- 
tings between the pairs r = 1,2 and r = 3,4 is much 
larger then splittings within the pairs. The latter is usu- 
ally not resolved experimentally. If we completely neglect 
second tunneling, we are left with only two states, but 
for consistency they are labeled as r = 1 and t = 3 states 

Gas. 

Low energy effective Hamiltonian has the form (atomic 
units are used throughout the paper): 

-ffcff = E T + H rot , (1) 

where E T is tunneling energy and H rot is Hamiltonian of 
the rigid asymmetric top, whose constants A, B, and C 
weakly depend on the quantum number r (see Table [I]) . 



TABLE I: Parameters of the effective Hamiltonian {TJ in GHz 
from Ref. fTl. 
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E T 


A T 


B T 


cv 


1 


0.0 


301.873 


26.193 


25.120 


3 


342.885 


301.585 


26.142 


25.201 



Dimensionless sensitivity coefficient of the transi 
tion uj to /x-variation is defined so that: 
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In order to find these coefficients for transitions described 
by Hamiltonian (JXJ) we need to know /i-dependence of 
the parameters from Table HI Rotational parameters A, 
B, and C scale linearly with /i. Thus, purely rotational 
transition has sensitivity Qp. ro t = I- Sensitivity of the 
tunneling energy E T to \x- variation can be estimated with 
the help of semiclassical Wentzel-Kramers-Brillouin ap- 
proximation. Following [19] . we can write: 

E T = 2 -^e~ s , (3) 

where S is the action over classically forbidden region and 
Eq is the zero point energy for the inversion mode. Ex- 
pression ([3]) gives the following sensitivity to /i-variation: 

Q = S±l + ^2 ; ( 4 ) 

2 2([/ max — Eq) 

where U max is the barrier hight. Numerical factor in the 
second fraction depends on the barrier shape. For the 
triangular and square barriers it is 2 times smaller and 
1.5 times larger respectively @. The factor i in Eq. (HJ 
corresponds to the parabolic barrier [9j . 

The potential for the tunneling coordinate 7 was found 
by Koput et al. [2pj in a form of Fourier expansion with 
U m ax = = 387 cm -1 . Numerical solution of the one 
dimensional Schrodinger equation for this potential gives 
E Q = 169 cm" 1 E T = 13.1 cnT 1 and Q^ T = 2.44. Using 
this value for Eq and tunneling energy from Table U we 
can find S from Eq. (|3|) and Q^.r from Eq. (gj: 

5 = 2.28, Q M , T =2.54. (5) 

If we scale potential from [20| to fit experimental tunnel- 
ing frequency from Table U the numerical solution gives 
Q^.t = 2.56. All three values for Q M , r agree to 5%. The 
semiclassical value corresponds to the experimental tun- 
neling frequency and is less sensitive to the details of the 
potential shape. Therefore, we use it in our calculations 
and assign it 5% uncertainty. 

We see, that tunneling energy is 2.5 times more sensi- 
tive to /i-variation than rotational energy. Sensitivity of 
the mixed tunneling- rotational transitions u) = E T ±Lo rot , 
is a weighted average of the tunneling and rotational con- 
tributions [21| : 

Qn = Qu.t i <3u,rot • (6) 

This sensitivity is further enhanced for the frequencies 
\u)\ < E T . 

Results of the numerical calculations with Hamiltonian 
([1]) are given in Table [TTJ In this table we list six transi- 
tions in the range 200 - 700 GHz, which were studied in 
Ref. [l6| (four of them were observed and other two were 
not) and transitions from the JPL database [22j] with fre- 
quencies below 100 GHz. 

Note that pure rotational transitions with At = for 
peroxide are not observed. Because of that all sensi- 
tivities in Table [TT] significantly deviate from unity. As 



TABLE II: Numerical calculation of the Q-factors for low fre- 
quency mixed transitions in peroxide using Hamiltonian (JT| 
and Eq. (01 . Experimental frequencies are taken from JPL 
Catalogue [221 ]. _E up is upper state energy in Kelvin. 



Jk a 


k c (r) 


-E up 




(MHz) 


Q 




upper 


lower 


(K) 


theory 


exper. 








Transitions below 100 GHz 






0o,o(3) 


-li,o(l) 


17 


14818.8 


14829.1 


+36 


.5(2.9) 


2i,i(l) 


- lo,i(3) 


21 


37537.0 


37518.28 


-13 


.0(1.2) 


!o,i(3) 




19 


67234.5 


67245.7 


+8 


.8(6) 


2o,2(3) 


-2i,a(l) 


24 


68365.3 


68385.0 


+8 


.7(6) 


3o,3(3J 


-3i l3 (l) 


31 


70057.4 


70090.2 


+8 


.5(6) 


An t 

^0,4 \o) 


- 4i, 4 (l) 


41 


72306.0 


72356.4 


+8 




Jo, 5 yd) 


-5i,s(l) 


53 


75104.6 


75177.4 


+8 


• U W 


6n fi(3) 


-6i.e(l) 


68 


78444.7 


78545.4 


+7 


7(6) 


3i,a(l) 


-2o, 2 (3) 


28 


90399.8 


90365.51 


-4 


.8(5) 




Transitions observed from ISM in Ref. [16] 




3o, 3 (3) 


-2i,i(l) 


31 


219163.2 


219166.9 


+3 


.4(2) 


61,5(1) 


-5o, 5 (3) 


66 


252063.6 


251914.68 


-1 


.1(2) 


4o,4(3) 


-3i, a (l) 


41 


268963.7 


268961.2 


+3 


.0(2) 


5o,s(3) 


-4i, 3 (l) 


53 


318237.7 


318222.5 


+2 


.7(1) 




Transitions attempted to observe in 


Ref. [16] 




5i, 4 (3) 


-6o, 8 (l) 


67 


318635.6 


318712.1 


+2 


.7(1) 


li,o(3) 


-0o,o(l) 


32 


670611.9 


670595.8 


+ 1 


.8(1) 



expected, lower frequency transitions have higher sensi- 
tivities. Transitions with the frequency \oj\ < E T fall 
in two categories. For transitions t = 3 — > r = 1 the 
tunneling energy is larger than rotational energy and 
= E T — w ro t. For such transitions Q-factors are posi- 
tive. For transitions r = f— >t = 3, \u>\ = ui I0 t — E T and 
Q-factors are negative. 

Let us discuss the accuracy of our estimates. Calcu- 
lated frequencies agree with experiment to 0.1%, or bet- 
ter in spite of the simplicity of Hamiltonian (JlJ. That 
means that centrifugal corrections to the rotational en- 
ergy are indeed unimportant for the transitions with low 
rotational quantum numbers considered here. Even in 
this case the scaling of the rotational energy with \i is 
not exactly linear as rotational constants of the effective 
Hamiltonian are averaged over the vibrational wave func- 
tions of the states r = 1,3. Respective corrections to the 
rotational sensitivity coefficients are of the order of 1% - 
2% [ll|] . Note that r-dependence of the rotational con- 
stants can be considered as centrifugal corrections to the 
tunneling energy 0. According to Refs. d |Tlj] the ac- 
curacy of the semiclassical expressions for tunneling (0 
[4]) is comparable, about 3%. Additional uncertainty is 
associated with the zero point energy Eq, which is not 
directly observable and we take it from calculation with 
potential [20j | . This potential was used in a number of 
later papers, including [23|, [24j. We estimate the accu- 
racy of the tunneling sensitivity coefficient ([5]) to be 5%, 
which agrees with numerical estimates made above. Un- 
certainties in rotational and tunneling sensitivities result 
in the final uncertainties for mixed transitions given in 
Table [TTJ These uncertainties are typically about 10%. 
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Such accuracy is sufficient for the analysis of the astro- 
physical spectra. 

To summarize, we have estimated sensitivity coeffi- 
cients for the low frequency mixed tunneling-rotational 
transitions in H2O2. Transitions below 100 GHz appear 
to be highly sensitive to /i-variation with sensitivity 
coefficients of both signs. Maximal relative sensitivity 
is found for transitions 14.8 and 37.5 GHz, where 
AQ M = +36.5 - (-13.0) « 50. This is comparable 
to the highest sensitivities in methanol [l(| [Til and 
an order of magnitude larger than for ammonia [ll l2l- 
The lines recently observed by Bergman et al. [16[ 
have frequencies above 200 GHz and significantly lower 
sensitivities. Three of the observed transitions have 
sensitivities, which are close to each other. However, 
the sensitivity of the fourth observed transition is 
significantly different. Transitions 219 GHz and 252 
GHz have relative sensitivity AQ M = 4.5. This is 



close to the sensitivity of the ammonia method, where 
= 3.5, but with the advantage that both lines 
belong to one species. This eliminates very important 
source of systematic errors caused by the difference 
in spatial and velocity distributions of species Q. 
We conclude, that microwave transitions in H2O2 
potentially can be used for /i-variation search in as- 
trophysics. Note, that high sensitivity transitions 
correspond to low rotational quantum numbers J. This 
makes peroxide a potential candidate for labora tory 
tests on /i-variation using molecular beam technique [25( ■ 
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